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Abstract 

Boundary S matrices for the boundary tricritical Ising field theory (TIM), both 
with and without supersymmetry, have previously been proposed. Here we provide 
support for these S matrices by showing that the corresponding boundary entropies 
are consistent with the expected boundary flows. We develop the fusion procedure for 
boundary RSOS models, with which we derive exact inversion identities for the TIM. 
We confirm the TBA description of nonsupersymmetric boundary flows of Lesage et 
al, and we obtain corresponding descriptions of supersymmetric boundary flows. 
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1 Introduction 



A well-known (but nevertheless, remarkable) feature of integrable quantum field theories in 
1 + 1 dimensions is that their exact bulk Q and boundary || scattering matrices can be 
found. However, such results are generally not obtained in a systematic way from the action; 
rather, one often relies on general principles (factorizability, unitarity, crossing, bootstrap, 
etc.) and educated guesses about symmetry, mass spectrum, etc. A case in point is the 
tricritical Ising field theory - i.e., the tricritical Ising conformal field theory (CFT) [^, [|. Q 
perturbed by the $(i i3 ) operator || . We shall refer to this field theory as the "tricritical Ising 
model " or TIM for short. The bulk S matrix was proposed in 0, and boundary S matrices 
were proposed in || [|. This field theory has several notable properties, which render it a 
very attractive toy model: it is unitary; it is supersymmetric; and it is one of the simplest 
examples of a model of massive kinks, whose scattering matrices are of RSOS [|K], [ll[] type. 



Moreover, the bulk and boundary soliton S matrices [12], [L3| of the N — 1 supersymmetric 
sine-Gordon model [14], |15| contain the corresponding TIM S matrices as one of the factors. 

A thermodynamic Bethe Ansatz (TBA) analysis [16] can provide a nontrivial check on a 



given bulk [T7L 18] or boundary [19, EO, EE scattering matrix. Indeed, the S matrices serve 



as the input of the "TBA machinery," whose output consists of certain data (central charge 
||. f|], boundary entropy |22|, ^3| ) which characterizes the corresponding CFT. For the TIM, 
a TBA check of the proposed bulk S matrix [7] was performed in ||18|| . 

One of the principal aims of this paper is to perform an analogous TBA check of the 
boundary S matrices which have been proposed in || ||. Such an analysis is technically 
nontrivial, since neither the bulk nor boundary S matrices are diagonal. As in the bulk 
case JT8f, the key step is the derivation of an exact inversion identity which is obeyed by an 



appropriate transfer matrix. For the boundary case considered here, the transfer matrix is 
of the "double-row" type p4 . 



A second aim of this paper is to develop the techniques for deriving the necessary inversion 
identity. We do this in an extended appendix, building on earlier work on fusion for vertex 



2g, |27| and RSOS ]28|, £9], |30[ models. The main idea is to formulate an RSOS open- 
chain fusion formula, and to show that the TIM fused transfer matrix is proportional to the 
identity matrix. 

A third aim of this paper is to derive TBA descriptions of TIM massless boundary flows. 
Let us recall p2, || H that the tricritical Ising CFT has a discrete set of (super) conformal 
boundary conditions. Boundary perturbations can lead to flows among these boundary 



conditions |q, 31, 32, B3L R41 RJ. A TBA description of the nonsupersymmetric flows was 



proposed in ]^T|] on the basis of an analogy with the Kondo problem. Here we give a 
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derivation of that TBA result, as well as the results for supersymmetric flows not considered 
in |yj , directly from the TIM scattering theory. (An alternative approach based on a lattice 
formulation of the TIM is considered in [Q. However, it seems that this approach cannot 
generate the boundary entropies.) 

The outline of this article is as follows. In Sec. 0, we review the bulk and boundary S ma- 
trices ]7|, [| which will serve as inputs for our TBA calculation. There are two boundary S 
matrices that are not supersymmetric; and there are two boundary S matrices which do have 
supersymmetry, which we call NS and R. We also briefly review the classification of (super) 
conformal boundary conditions, certain pairs of which are connected by boundary flows. In 
Sec. |3|, we carry out the first step of the TBA program, which consists of constructing the 
so-called Yang matrix [[35|] and relating it to a commuting transfer matrix. For the problem 
at hand, we require a boundary RSOS version of the Yang matrix, which is an interesting 
generalization of the known case of periodic boundary conditions. In Sec. |], we use an exact 
inversion identity to determine the eigenvalues of the transfer matrix in terms of roots of 
certain Bethe Ansatz equations. We restrict our attention here to the NS case. In Sec. |5] 
we use these results to derive the TBA equations and boundary entropy. Moreover, we find 
massless scaling limits which correspond to boundary flows, both for the NS case and the 
nonsupersymmetric cases. In Sec. ^] we briefly discuss the R case, which is closely related 
(in fact, dual) to the NS case. Our conclusions are presented in Sec. [7|. In an Appendix, we 
give a brief account of the fusion procedure for RSOS models with boundary, and provide 
the derivation of the TIM inversion identity. 



2 TIM scattering theory 

We briefly review in this Section some pertinent results on the TIM scattering theory. We 
first define the bulk model as a perturbed bulk CFT, and give the bulk S matrix fTfl. We 
then enumerate the possible (super) conformal boundary conditions, and give the boundary 
S matrices which have been proposed 0, || to describe certain perturbations of some of these 
boundary conditions. Two of the boundary S matrices do not have supersymmetry, and two 
of them do. Many of the notations used in this paper are introduced in this Section. 

2.1 Bulk 

The bulk TIM is defined by the "action" |7j 

/*oo /*oo 

A = A MW5) + \ dy dx^^z^x.y), A<0, (2.1) 
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Figure 1: Bulk S matrix S c a £(0i - 9 2 ) 



where AjK( 4( / 5 ) is the action for the tricritical Ising CFT (i.e., the minimal unitary model 
A4(4/5) with central charge c = jq), and $^3 3-, is the spinless (1 ,3) primary field of this 
CFT with dimensions (| , |). Moreover, A is a bulk parameter with dimension length - ^. We 
restrict our attention to the case A < 0, for which there is a three-fold vacuum degeneracy, 
and the spectrum consists of massive (mass m > 0) kinks K a ^{9) that separate neighboring 
vacua, a , b G { — 1 ,0,1} with \a — b\ = 1. Multi-kink states 



K, 



K. 



a.2 ,b2\y2, 



must obey the adjacency conditions b\ = a 2 , etc. 

The two-kink S matrix S° f(9) is defined by the relation (see Figure H) 

K atC (0i) K c>b (6 2 ) = Y, S a " ^) K a>d (B 9 ) Kd 



(2.2) 



The nonzero matrix elements are given by J7|, [L8| [] 



syAo) 



°(e)sy a ,{e) 



(2.3) 



where a ,a' G { — 1 ,+1}, 7 = n- In 2, and the "reduced" matrix elements f{6) are given 
by 



cosh • 



S° 7(0) = -isinh 







1 



cosh -(9 — in) 
4 



S° a (6) = V2cosh±(6 + zn) 



Finally, a (9) is a function which obeys 

a(6) = - 9) , a{9) a{-9) 



cosh : 



(2.4) 



(2.5) 



x It is noted in j7j that this S matrix is essentially the solution of the star-triangle equation corresponding 
to the critical Ising lattice model [|o[ |ll) . 
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and has no poles in the physical strip < Sm 9 < n. A useful integral representation for 
this function is 

(e) _ (, f°°dt sin(g)sinh(f) 

a[ ' V2 sinh(±(# - in)) \ Jo t sinh(2t) cosh(|) 

This is a "reduction" of the well-known integral representation for the factor U{6) of the 
sine-Gordon S matrix |l| with ^ = ~. 

Zamolodchikov has shown in |7[] that this S 1 matrix "commutes" with supersymmetry 
charges Q and Q, which obey the N — 1 supersymmetry algebra with topological charge. 
The S matrix also commutes with the spin- reversal operator T, which is defined by 

T K aim {9l) Ka 2 ,a 3 (0 2 ) . . ■ K aN jaN+1 (8 N ) 

= K. ai ,_ a2 (9 1 ) K. a2 ,_ a3 (0 2 ) . . . K_ aiV ,_ aiV+1 (^) . (2.7) 
Further properties of the S matrix are listed in Appendix [A]. 




2.2 Boundary 

Although the three vacua —1,0,-1-1 are degenerate in the bulk, these vacua do not necessar- 
ily remain degenerate at the boundary. Chim || has identified the six conformal boundary 
conditions (CBC) |22] of the tricritical Ising CFT as follows: for the boundary conditions 
(— ) , (0) , (+), the order parameter is fixed at the boundary to the vacua —1,0,-1-1, re- 
spectively. For the boundary condition (—0), the vacua —1 and are degenerate at the 
boundary; hence, the order parameter at the boundary may be in either of these two vacua. 
Similarly, for the boundary condition (0+), the and +1 vacua are degenerate at the bound- 
ary. Finally, for the boundary condition (d), all three vacua —1,0, +1 are degenerate at the 
boundary (as well as in the bulk); i.e., the order parameter at the boundary may be in any 
of the three vacua. The corresponding g factors |23j are given by [|J 



g {d) = V2ifC , # { _o) = #(0+) = rfC , 

#(0) = V2C , # ( _) = g {+) = C , (2.1 



where 



, sin f / sin ^ , 

It is argued in that the conformal boundary conditions (—)&(+), (— 0)&(0+), (0) and 
(d) are in fact superconformal. Notice that the first two of these super conformal boundary 
conditions correspond to superpositions of "pure" Cardy states. 
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Figure 2: Boundary S matrix Ra [0). 



b 



We shall consider separately integrable perturbations of both conformal and superconfor- 
mal boundary conditions, resulting in models without and with supersymmetry, respectively. 
We assume || that also in the perturbed theory the boundary can have (at most) three pos- 
sible states, corresponding to the three different vacua, which are created by the boundary 
operator B a with a £ {—1 ,0,1}. Multi-kink states have the form 



K ai , a2 (0i) K a2ta3 (9 2 )...K aN!a (9 N ) B, 



The kink boundary S matrix R a (#) is defined by the relation (see Figure 0) 



b 




(2.10) 



c 



2.2.1 Non-supersymmetric cases 



Chim |]] has considered the TIM on the half-line x < corresponding to an integrable 
perturbation of the CBC (—0). The model is defined by the action 



.4 





A < . 



(2.11) 



The last term is the boundary perturbation. It involves the boundary primary field <j>^ ^_ ^ 
with dimension A( 13 ) = | which acts on the CBC (—0). | Moreover, h is a boundary 

2 

parameter which has dimensions length - 5. 

The boundary S matrix which has been proposed || for this model has the following 
nonzero matrix elements 

R°l[ (o,0 = p(0,o£o~* 

R«° Q (0,0 = M(e,Z)R*° Q (0,0, (2.12) 

_ c 

where the reduced matrix elements R a (0,0 are given by 



i> 



-1 



R*^ (0,0 = 1, 

R ±1^(6,0 = cos|±isinh-. (2.13) 

The parameter £ is related in some way to the boundary parameter h appearing in the action 
( [2. 11| ). The function P(8 ,0 is given by 

P(9,O = P CDD (0,O P min (0), (2.14) 
where P CDD (6 ,0 is the CDD factor 

peon si^-isinhe 
sin 4 + 1 sinn 9 

which has a pole at 9 = iC,, and P m in(0) is the minimal solution of the equations 

Prmn(0)P m n(-0) = 1 , P min (y - 0) = V^ 2 ^ COsh(^ - j)(T (29) P min (j + 0) , (2.16) 

with no poles in the physical strip < Sm 9 < |. We find that it has the integral represen- 
tation 

Finally, the function M{9 , £) is given by 

M{9,Q = e 2 ^ a(9 - i£) a(9 + i£) P(0 , f) . (2.18) 

2 In general, boundary operators (f) a and which act on conformal boundary conditions a and b commute; 
i.e., their operator product expansion with each other is zero. Such operators have recently been studied in 
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There is a similar model corresponding to a perturbation of the CBC (0+). The boundary 



-i 



S matrix for this case is the same as the one given above, except that R o (9 , £) = 0, and 



Rv (e,O = P(0,Z)R\ (9,0 



(2.19) 



with R o (6 ( £) = 1. Neither of these two models has supersymmetry. 



2.2.2 Supersymmetric cases 

Super symmetric perturbations of the tricritical Ising boundary CFT with two different su- 
perconformal boundary conditions (namely, (— 0)&(0+) and (d)) are considered in ||. We 
refer to these two NS and R, respectively, since these are the sectors to which the 

corresponding boundary states belong. 

NS case 

The NS case corresponds to a perturbation of the boundary condition (— 0)&(0+), with 
action 

A = Am(4/5)+(-0)&(0+) 



/oo (-0 
dy / dx $/3 i |)(x , y) 
-oo J —oo 



hi dy U 

J —oo 



(l),(-o) 



(y) 



'(|),(0+)(2/) ) > A < °- 



(2.20) 



The proposed boundary S matrix is the "direct sum" of the boundary S matrices given in 



Sec. 2.2.1 for the perturbations of (—0) and (0+). That is, the nonzero matrix elements are 
given by 



Ro (9,0 = P(9,0R° (0,0, 

a a 

R*° (9,0 = M(9,0 R'l (0,0 



(2.21) 



where 



Ro (9,0 = 1 



o £ 9 

R±i (9,0 = cos - ± % sinh - 
o v ' s; 2 2 



(2.22) 



and the functions P(9,0 an d M(9 ,0 are given by Eqs. (|2.14|) and ( p,18|) , respectively. 
This boundary S matrix "commutes" with the supersymmetry charge 



Q = Q + Q + 2cos(^) v / ^r, 



(2.23) 
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where T is the spin- reversal operator (|2.7| ). 



R case 

For the R case, which corresponds to a perturbation of the boundary condition (d), the 
action is given by the image of ( [2.20|) under duality transformation. The proposed boundary 
S matrix has the following nonzero matrix elements 

Ro a (9,0 = N(9^)R° a (0,0, 

a a 

R^ a (o,0 = N(6,0R°~ a (0,0, 

a a 

R*° q (6,0 = R(9,0R°l (0,0, (2-24) 

c 

where the reduced matrix elements R a (6,0 are given by 



i, 



Ro a (9,0 = cos| 



Ro^(6,0 = -*rsinh^, R o ^ (6 ,0 = -- sinh °- , 

R°° o (0,0 = 1, (2-25) 

and r is a parameter which presumably is related in some way to the boundary parameter 
h, as is £. Moreover, the functions N(6 ,0 and R(9 ,0 are given by 

N(6,0 = e- 2 ^ d M(6,0, 

R(9,0 = e 2 ^ e P(6,0- (2.26) 
This boundary S matrix "commutes" with the supersymmetry charge 



— Ait £ 

Q = Q-Q + - -cos(i)V^r. (2.27) 

1 — r z 2 



In contrast to the NS case, here the matrix R a (8,0 does not vanish for b ^ c; i.e., it is 

b 

not "diagonal." 

The parameter r can be set to unity by an appropriate gauge transformation [[| of the 
kink operators, which corresponds to adding a total derivative term to the boundary action 
that restores spin-reversal symmetry. This limiting case, for which the supersymmetry charge 
02.271) reduces to T, was considered earlier in ||. 
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3 Yang matrix and transfer matrix 



The first step of the TBA program is to formulate the "Yang matrix" [35 and relate it to an 
appropriate commuting transfer matrix. Since it is not obvious how to do this for the case 
of boundaries, we begin by reviewing the case [ITS] of periodic boundary conditions. f\ 



3.1 Closed-chain transfer matrix 



Following flTg| , we consider N kinks of mass m with real rapidities 9\ , . . . , 9^ and two-kink 
S matrix S% f(9) in a periodic box of length L >> — . We impose the periodicity condition 

e iLM ^K aua2 {e x )K a2m {6 2 ) . ..K aN _^ aN {9 N ^)K aNAl {9 N ) 

= K a2 , az {9 2 ) . ..K aN>ai {9 N )K ai>a2 {9i) . (3.1) 

Commuting the kink operator K ai ia2 (0i) on the LHS past the others using the relation (|2.2| ), 
we obtain 

elLmM J- {S^ J(0 a - 9 2 )S% Jft -9 3 )... S%zl Z"^ - °n-i)SZ_X(0i - On) 

d'2 ,rfjv 

K ai id2 (9 2 )K d2 id3 (9 3 )... K aN _i , dN (9 N )K dN , ai (6 1 ) } = K a2 , a3 (9 2 ) . . . K aN >ai (9 N )K ai , 02 (9 1 ) . 

(3.2) 

Multiplying both sides by the "wavefunction" \l/ ai ■•• ctiV , summing over a\ , . . . , ajv, and rela- 
beling indices appropriately, we obtain the Yang equation for kink 1 



giLm sinh 9\ 



Yl a a\-X ^ a,1 '" a ' N = y ai ~ aN , (3.3) 



where Ym is the Yang matrix 

y w zz = 5 < s i Z& - e * )s i Z 9 * s t-i - 9n - i)s «* ^ - 9n) ■ (3 - 4) 

There are similar equations, and corresponding matrices Yha, for the other kinks k = 
2,3,... ,N. 

The objective is to diagonalize Yn^. The key to this problem is to relate Yha to an inho- 
mogeneous closed-chain transfer matrix, for which there are well-developed diagonalization 
techniques. To this end, we consider the transfer matrix 



3 The analysis presented here for RSOS-type S matrices is parallel to the one given in pl| for vertex-type 
S matrices. 
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with inhomogeneities 9\ , . . . ,9^. Because the S matrix satisfies the Yang-Baxter equation 
( |A.5| ), the transfer matrix commutes for different values of 9 



... ,9 n ),t{9%,... ,9 n )} = 0. (3.6) 

Let us now evaluate this transfer matrix at 9 = 9\. Using the fact that the S matrix at 
zero rapidity is given by (|A.6|) , we immediately obtain t(9i\9i , . . . , 9n) = 111 general, 
we have 

Y (k) = r(9 k \9 1} ... ,9 N ), k = l,... ,N. (3.7) 

This is the sought-after relation. In order to diagonalize the Yang matrices F(m, it suffices to 
diagonalize the commuting closed-chain transfer matrix t{9\9\ ,.. . ,9^). That calculation, 



as well as the corresponding bulk TBA analysis, is described for the TIM in [IS 



3.2 Open-chain transfer matrix 

We turn now to the case with boundaries, which is our primary interest here. We therefore 
consider N kinks of mass m with real, positive rapidities 9\ , . . . , 9^ in an interval of length 

c 

L » — , with bulk S matrix S° f(9) and boundary S matrix R a [9 ,£). In analogy with 
(^.1[), we propose the formal relation 

e 2iLm S inM 1B + ^ , a3 (9 2 ) . . . K att _, ,a N (9 N ^)K aN , ajv+1 (M 

= B X K ai,a2{9i)K a2ia3 {9 2 ) ■ ■ ■K aN _ 1 1 a N (9 N ^ 1 )K aN >aN+1 (9 N ) B~ n+i , (3.8) 

where now there are two boundary operators corresponding to the left and right bound- 
aries, with (cf, Eq. ( gTop ) Q 



K a , b (9)B b = J2 R \( 9 ^-) K "A-0)B;, (3.9) 

c 

B+K b , a (9) = ^2 B c K c>a (-9) Ra° (-9, £+) . (3.10) 



4 Our convention for matrix multiplication is given by 

V )a[...a' N Z_^i a[...a' N a'{...a'^ ' 



5 The relations ( |3.9| ) and (3.10) are consistent in that both lead to the same boundary Yang-Baxter 
equation ( |A.1Q ). 
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Note that for each boundary operator there is a corresponding boundary parameter £±. 
By moving the kink operator with rapidity 9\ on the LHS of (|3.8| ) to the far right using (|2.2|) , 



reflecting it from the right boundary using ( |3.9| ), moving it to the far left using again ( |2.2| ). 
and finally reflecting it from the left boundary using ( |3.10|) , we arrive at the Yang equation 
for kink 1 



JliLm sinh 



h \ ~" y ai—a-N+i \ua[... a' N+1 _ ^ai—ajv+i 
(1) a[...a' N+1 



(3.11) 



JV+l 



where the Yang matrix Km is given by 



Y ai...a N+ i 
(1) a[...a / N+1 



X 



i? , £_) Sg,, S +1 (fi + 9 N )... 5* J + 2 ) ) . (3.12) 



1^2 ,<ijv 
a JV+l 
J JV+1 

There are similar matrices Yqa for the other kinks. In analogy with the case of periodic 
boundary conditions, the key to diagonalizing the Yang matrix is to relate it to an inhomo- 
geneous open-chain transfer matrix 



fj~y +1 (6\6i 



e 



i 1 ...u JV + 1 



N) 



[in 



S C 4 



C a J i 



"jv "JV+l 



Nj 



X i?a 



JV+l 



c^ 1 (e, e_) ^ ajv (^ + . . . s$ r + 

i 1 1 1 



^JV+l 



which commutes for different values of 9 

\t(9\9\ , • • • , 0/v) , t{9'\9i , . . . , #tv) 



0. 



(3.13) 



(3.14) 



The transfer matrix (|3.13|) is an RSOS version |T^, |3(J] of the Sklyanin p4] vertex-type 
transfer matrix. Using the relations (|A.6|) , (|A.9|) and (|A.2 ), one can show that 



K 



(fe) 







A: 



(3.15) 



Hence, in order to diagonalize the Yang matrices Ko.), it suffices to diagonalize the open-chain 
transfer matrix t{9\9\ , . . . , 9^). Indeed, let \l/(#i , . . . , 9^) be an eigenvector of the transfer 
matrix with corresponding eigenvalue A(9\9i , . . . ,9jy), 



t(9\9 u ... ,9 N ) ,9 N ) = A(9\9 U ... ,9 N ) (0 X , . . . , 



A'/ 



(3.16) 



The eigenvector is independent of 9 by virtue of the commutativity property ( |3.14j ). With 
the help of the result (|3.15|) , the Yang equation (|3.11|) implies 



3 2iLm sinh 8^ 



A(9 k \9i , . . . , 9 N ) 



k = 1 N . 



(3.17) 
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4 Inversion identity and transfer-matrix eigenvalues: 
NS case 



We turn now to the problem of determining the eigenvalues of the inhomogeneous open- 



chain transfer matrix (|3.13|) . As for the closed chain [ IS ] , our approach is to derive an exact 
inversion identity. For definiteness, we treat here the NS case. (See Sec. |2.2.2| .) The results 
for the R case, which are closely related to those for the NS case, are presented in Sec. |6|. 

Instead of working with the full ("dressed") transfer matrix ( |3.13| ), it is convenient (see 
Footnote |^ below) to work instead with the reduced ( "bare" ) transfer matrix t, which is 
constructed from the reduced bulk and boundary S matrices, 



a-, ...a 



i—^n+i 



x Ra' k (9,^)S a ^ N+ y + 9 N )...S:r a l(9 + 9 1 )}. (4.1) 

a N+l J 

It is also convenient to define the following four "sectors" : 

N = even — Sector I : a x ,a[ , a N+ i , a' N+1 G { — 1 , +1} , 

N = even — Sector 77 : a\ = a[ = a^ + i = a' N+1 = , 

N= odd — Sector J : a\ , a\ G { — 1 , +1} , cln+i = a N+i = > 

N= odd — Sector// : ai = = , a^+i ,a' N+1 E { — 1 ,+1} . (4.2) 

The nonzero matrix elements of the transfer matrix lie exclusively in these sectors. For a 
given parity of N (i.e., even or odd), the transfer matrix decomposes into two blocks along 



the diagonal corresponding to sectors I and II. For the NS case (|2.21|), (|2.22|), the relation 



between the full transfer matrix and the reduced transfer matrix is given by 

t^{e\e x ,...,9 N ) = w {a \9) ^(919, ,...,e N ), (4.3) 

where a runs over the four sectors ( |4.2|) , and w^ a \9) is given by 



N 

w {a \9) = Y[a(9 - 9 3 )a{9 + 63) x < 
3=1 



{ P(z7r-0,£ + ) 

M(m-e,z + ) M(e,t-) 

e- 2h > e P(in-9,£ + ) M(0,£_) 
[ e 2 ^ e M{iix-9^ + ) F(0,£_) 



(4.4) 
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respectively. The latter can be brought to the form 



N 



a(26) cosh f 1 = 1 



\\a{9 -9 j )a{9 + 9^ x < 



{ e^ 6 P(9,C+) P(e,£J) 
\e^ e M{9^ + ) M(9,C- 
P{9,Z+) M(0,£_) 
{ ±M(0,£ + ) P(0,£_) 



(4.5) 



with the help of the crossing properties ( [2.5| ), ( |2.16| ). 

Using the fusion procedure, we show in Appendix [B] that the reduced transfer matrix 
obeys the inversion identity 



# a \9\9 x ,...,9 N ) i^\9 + i7r|0! ,...,9 N ) = f^(9) I™ 
where a runs over the four sectors ([4.2[), and f^ a '(9) is given by 



(4.6) 



f {a) {0) 



1 



cosh 9 



fi a \9) cosh 2 u - J] cosh(^ - 9 3 )) cosh(^ + B,)) 



9 



N 



1 



where 



+ /i a) (0) sinh 2 t f[ sinh(i(0 - %)) sinh(i(0 + %)) 

1 

(cosh^ ± cos £_) (cosh 6* ± cos£_|_) 
| (cosh ± cos£_) 



(4.7) 



f { ±\0) = { 



2 (cosh 9 ± cos£+) 

respectively. This inversion identity is one of the main results of this paper. We have checked 
it explicitly up to N = 4. 

In addition to the inversion identity, we can establish certain further properties of the 
transfer matrix which are needed to determine its eigenvalues. Namely, periodicity [] 



t{9 + 2ztt|0i , . . . ) 9 N ) = t{9\9 l) ... ,9 



A' J 



crossing 

t(i7r-0|0 ls ... ,9 N ) =t(0|0!,... ,9 N ), 

and asymptotic behavior for large 9 

^ a \9\9 u ... ,9 N ) ~ z (a) {9) for 0^oo 

6 This is not the case for the full transfer matrix t(9\9\ , . . . , On). 



(4.9) 



(4.10) 



(4.11) 
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where a runs over the four sectors (f4.2|) , and z^ a > is given by 



z (a) (9) = { 



N_ 
\ 2 



(^ai .ajv+i <^ai 



-ajv+i , 



'JV+l 
2 



+ 1 



(4.12) 



-1 - ^ai ,1> 



N+l 

2[-^-| (5 ajv+1 -i - 5 aiV+1 ,i) 



respectively. 

Acting with the above relations on an eigenvector ^{9% , . . . , 9^) of the (reduced) transfer 



matrix 



i{e\e u ... ,e N ) ^ , . . . ,e N ) = K{e\e x , . . . ,e N ) *{9 U ...,9 



(4.13) 



we obtain corresponding relations for the eigenvalues A^ a \9\9i , . . . , 6n) in the various sec- 
tors, 

A^\9\9 U ... ,9 N ) A^(9 + i7r\9 1 ,... ,9 N ) = f^(9), (4.14) 
A<«>(0 + 2i7r|0i,... ,9 N ) = A^ a \9\9 1 ,...,9 N ), (4.15) 
K^(i>K-9\9 1 ,... ,9 N ) = AW^,...,^), (4.16) 
A^(0|0i,... ,Sjv) ~ z {a \9) for tf^oo. (4.17) 



The periodicity, crossing and asymptotic behavior requirements of the eigenvalues ( |4.15|) 
( |4.17| ) are fulfilled by the Ansatz 



A^^ , . . . , 9 N ) = c< a > sinh(^(0 - Uj )) cosh(i(0 + Uj )) , 



(4.18) 



where and are given by 



r ±i 

2 

±1 
±2 



d {a) 



( N 

2 

N 



. +1 

JV+l 
2 

JV+l 



(4.19) 



< 2 



respectively. The parameters {uj} appearing in the Ansatz ( |4 . 1 8| ) are evidently roots of the 
eigenvalues, A^ a \uj\9i , . . . , On) = 0. It follows from the inversion identity ([4. 14|) that {uj} 
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are also roots of the function f^ a \9), i.e., f( a '(uj) = 0. We conclude from (O) that {uj} 
are solutions of the set of equations 



f (uj) sinh2 f TT sinh(|( Wj - 9 k )) sinh(|(^ + 6 k )) _ {a) 

/W( Ui ) cosh 2 f W cosh(|( M , - 0,)) cosh(|( u , + e k )) ' J ' • • • ' ' 1 • UJ 

to which we refer as "Bethe Ansatz" equations. 

The periodicity property (|4.15|) implies that we can restrict the roots Uj of A^ a '(9) to the 
interval 



— 7T < 3m Uj < 7T . 



(4.21) 



We now demonstrate that all the roots Uj have the form Xj ± y with Xj real. Indeed, we 
observe that f^ a \9) has the properties [] 



[f( a \9)]* = f {a \d* 



(«), 



(4.22) 



where * denotes complex conjugation. These two properties imply that if Uj is a root of 
f {a) (9), then so are u* and Uj =F in, respectively. Since u* ^ Uj, then u* = Uj =F m. Hence, 



3m Uj = ±y- 



In view of the above, we set 



with real and €j = ±1. The eigenvalues ( |4.18| ) are specified by {xj ,€j}, j = 1 



(4.23) 



similarly to the bulk case [jrj|. Hence, we can rewrite the expression for the eigenvalues as 



, . . . , e N ) = c w J] " A l 2) (^ + 



(4.24) 



where 



?7T 



A«(0) = sinh 1(0 --ej 



A( 2 )(^) = -zcoshi(^ + ^ e ) 



(4.25) 



Moreover, we can rewrite the Bethe Ansatz equations ( |4.20| ) in terms of Xj (they do not 
depend on 6j) as 



sinh 



27T \ N 



n 



sinh( J 2 fc — y) sinh( 



2 4 - 



sin h(^ + |) sinh (£2±^ + 1; 



1, 



J = l,...,^) 



(4.26) 



7 We assume here that and £± are real. 
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where 



and 



sinh x — i cos £ 
sinh x + i cos £ 



(4.27) 



(4.28) 



(As always, a runs over the four sectors Q4.2|) .) Notice that ( 4.26Q is invariant under 
—Xj. Moreover, following J37, [38], we assume that the root Xj = corresponds to an 
eigenvector with zero norm. Hence, we restrict to solutions with Xj > 0. 

To summarize this Section, the eigenvalues ( |3.16| ) of the full transfer matrix for the NS 
case of the TIM are given by 



A (a) (#l#i,--- ,9 N ) = w {a \9) A (a) (^|^ 







N) 



(4.29) 



where w^ a \6) is given by (|4.5|), A^ aS> {6\6\ , . . . ,9n) is given by ( |4.24| ), and {xj} are positive 
solutions of the Bethe Ansatz equations (|4.26|). 



5 TBA analysis 

Having obtained the eigenvalues of the transfer matrix and the Bethe Ansatz equations, we 
can proceed to the derivation of the corresponding TBA equations and boundary entropy. 



Following B, 19] we consider the partition function Z H of the system on a cylinder of length 



L and circumference R with left/right boundary conditions denoted by ±. It is given by 

Z + _ = tie- RB +-=e- RF 
= {B + \e~ LHp \B_} 

w (B + \0}(0\B_}e- LEo forL^oo. (5.1) 

In the first line, Euclidean time evolves along the circumference of the cylinder, and H + _ 
is the Hamiltonian for the system with spatial boundary conditions ±. In the second line, 
time evolves parallel to the axis of the cylinder, Hp is the Hamiltonian for the system 
with periodic boundary conditions, and \B±) are boundary states which encode initial/final 
(temporal) conditions. In the third line, we consider the limit L — > oo; the state |0) is the 
ground state of Hp, and E is the corresponding eigenvalue. The quantity ln(5 + |0)(0|S_) 
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is the sought-after boundary entropy [23, 19]. Taking the logarithm of the above expressions 
for the partition function, one obtains 



-RF &-LE + ki(B + \0)(Q\B-) . (5.2) 

Whereas the free energy F has a leading contribution which is of order L, we seek here the 
subleading correction which is of order 1. 



5.1 Thermodynamic limit 

We proceed to compute F using the TBA approach |l6|-[pTl. To this end, we introduce the 



densities P±{9) of "magnons" , i.e., of real Bethe Ansatz roots {xj} with ej = ±1, respectively; 
and also the densities pi{6) and p{9) of particles {9k} and holes, respectively. Computing 
the imaginary part of the logarithmic derivative of the "magnonic" Bethe Ansatz equations 
O), we obtain Q 

1 f°° 

p + {6) + p_{6) = — de' Pl {e') me-e') + <s>(e + e')} 

^ Jo 



where 



00 \smh(i + f )J coshfl' 

= ^aml. 0(».t)= 4C :t C ° 3 |l - (5.4) 
ad cos 24 + cosh 28 

In the final equality, we have used the expression ( |4.28|) for Q(9,£), and we have assumed 
that £ is real. We present here the results for N = even - Sector II, from which the results 
for the other sectors (see ( |4.27| )) can be read off by inspection. Defining pi(6) for negative 
values of 9 to be equal to pi(|0|), we obtain the final form 

P+(0) + P-{6) = -L ( Pl * $) (9) + ^ [m + %(0) + Ve_(9)] , (5.5) 



where * denotes convolution 



/oo 
iff f(9 - 9')g{9') . (5.6) 
-00 



8 There is a contribution — ^^r^W which originates from the exclusion |37], |38[ of the Bethe Ansatz root 



Xa = 0. 
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Computing the imaginary part of the logarithmic derivative of the Yang equation ( p. 17 ) 
using the result (|4.29 ) for the eigenvalue, we obtain (again for N = even - Sector II) 

Pl {6) + p{9) = ^-{2mcosh0 + ^ dff p x {9') [$ CT (0 - 9') + $ CT (0 + 0')] 



+ 



dff 



P + (6') u [ l\e - 6') + § 



(2), 



(2), 



+ 



1 

L 



-2 7 - $,(0) - 2^(20) + -9m lnM(0,f + ) + — 3m lnM(0,£_; 

at) ou 



(5.7) 



where 



d 1 /v 

— 3m In a{9) = - 



dk 



cosh 2 if 2vr sinh 9 ' 



*±(0) = ^m InA^(^), / = 1,2 



5 



,(0/ 



(5i 



and A± (0) are introduced in ( |4.25|) . Using the facts $± (0) = $± (0) = ±|<&(0), and defining 
P±{9) for negative values of 6* to be equal to P±(|0|), we obtain 

777 I ! 

Pl (9) + p(0) = _cosh0 + -(p 1 *$ (J )(0) + -((P + -P_)*$)(0) 

1 r a a 

+ — -2 1 -$ a (9)-2$ a (29) + — In M(9,t+) + —Zm In M(9,^)\5.9) 
ZttL I od o9 

Using ( |5.5| ) to eliminate P_, and ( |2.14| ), ( |2.18| ) to separate the various factors in M{6 , £), we 
obtain 



rn 



7l 



cosh 



1 

2^ 



1 



P + * $ H p x * $ CT — — $ * $ 



2tt 



1 



4tt 



2ttL 



9 



2 -3m lnP mm (0) - $ CT (20) - $ ff (0) + 7 + $ CT - — $ * $ 



<90 



4tt 



+ $ 



47T 



+ lnP OTD (0,£ + ) + lnP c ^( 



(5.10) 



Using the identity 



$ CT (0) - — ($*$) (0) = O. 

47T 



(5.11) 



as well as the identities 



$ IT (e-iZ) + $ (T (e + iZ)- — (Vt*$)(9) = o, 

lnP min (0)-$ CT (20)-$ a (0)+ 7 = -U>{9) 



we arrive at the final simple result 



pi(0) + = - coshfl + — (P + * $) (0) + 

7T Z7T Z7T.L 



where 



^ rx. r^rnn, n ^ 4 sin £ cosh 

= ™3m lnP cm? 0,£ = ^ i^TT 

<70 cos 2£ — cosh 20 



(5.12) 



(5.13) 



(5.14) 



The result ( |5.13| ) holds in fact for all four sectors. 

The thermodynamic limit of the magnonic Bethe Ansatz equations and the Yang equa- 
tions, given by (|5.5|) and (|5.13|) , respectively, are the main results of this subsection. 



5.2 TBA equations and boundary entropy 

The free energy F is given by 

F = E-TS, 
where the temperature is T = the energy E is 

L rco 

m cosh 6k = 

fc=i 



A' 



d9 pi(6)mcosh8 . 



(5.15) 



(5.16) 



E = m cosh 9k 

and the entropy S is |TB|, |T8| 

S = — J d9 |(pi + p) ln(pi + p) - pilnpi - plnp 

+ (P+ + P_) ln(P+ + P_) - P+ In P+ - P In P } . 
Extremizing the free energy (5F = 0) subject to the constraints 

5P- = -5P + + ^-5p! * $ , <5p = -8px + t^5P+ * $ , 

Z7T 27T 

(which follow from Eqs. ( |5.5|) , ( |5.13| ), respectively) we obtain a set of TBA equations which 



(5.17) 



(5.1* 



is the same as for the case of periodic boundary conditions [[Lqj 



tcosh0 = e x (0) + — ($*L 2 )(0), 
= e 2 (0) + -^($*L 1 )(0), 

Z7T 



(5.19) 
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where 



U{9) 



In (l + e 

Pi 



e 2 = In 



mR . 
Pi 



(5.20) 



We next evaluate F using also the constraints Q5.5[) , ( |5.13|) and the TBA equations. From 
the boundary (order 1) contribution, we obtain (see Eq. (|5.2j)) the boundary entropy 



ln(B + |0><0|fi_) 



1 

47T 



tZ0 



+ [f) + %W + $ ( _(S)]l 2 (e) . 



(5.22) 



In particular, the dependence of the boundary entropy of a single boundary on the boundary 
parameter £ is given by 



(5.23) 



where the kernels k§(#) and ^(0) are given in Eqs. ( 5.14]) and ( |5.4| ), respectively. This 
expression for the boundary entropy for the NS case of the TIM is another of the main 
results of this paper. 



5.3 Massless boundary flows 

We now consider the case of massless boundary flow. That is, we consider the bulk massless 
scaling limit 

m = /in , 9 = # =f In - , n -> , (5.24) 

where /i and 9 are finite, which implies E = fie ±e , P = ±fie ±0 . There are two nontrivial 
scaling limits of the boundary parameter, which we consider in turn. As we shall see, these 
two limits correspond to the boundary flows (-0)&(0+) -> (-)&(+) and (-0)&(0+) -> 2(0), 
respectively. 

9 Taking into account all the sectors, the last term in ( [5.22| ) should be replaced by 
[$(0) + VlK a )(0,£±)] L 2 {6), where a runs over the four sectors (fi~^) , and V?( a )(0,£±) is given by 

:;: ( ; ) + * s - (, ' , , (5-2i) 



*£+(«) 



respectively. 
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5.3.1 The boundary flow (-0)&(0+) 



-)&(+) 



Let us first consider the scaling limit 



n 







(5.25) 



7T 

\-l{. J , 

2 y 2 

where the boundary scale 9 b is finite. For the sign — in the limit ( 5.24)) , the CDD factor has 
a nontrivial limit 



P 



CDD, 



i sinh(- 



4 > 



sinh(- 



and therefore, so does the corresponding kernel ( |5.14| ) 

*e(9) -> $((9 - B ) . 



(5.26) 



(5.27) 



On the other hand, the factor Q(9,£) ( |4.28 ) becomes real in this limit; hence, the corre- 
sponding kernel ^$(9) (|5.4|) vanishes. The result ( |5.23|) for the boundary entropy therefore 
implies 



In q = — 



d§ §{9-9 B ) U{9) 



(5.28) 



where 6j(0) = €i{9 - lnf), and L^O) = ln(l + e~^ e) ). The factor of 2 appearing in fl5.28|) 
accounts for the contribution from the sign + in the limit ( 5.24j ), corresponding to the fact 
that right-movers and left-movers give equal contributions to the boundary entropy. In the 
UV limit 9b — » — oo, the integrand is nonvanishing for 9 — > — oo; similarly, the IR limit 
#b — ► oo requires ^ — > oo. Using the results Li(— oo) = In (|(3 + \/5)), Li(oo) = which 
follow from the TBA Eqs. ( 5.19 ), we conclude from (|5.28|) that 

MV 



g 



g 



IR 



This is precisely the ratio of g factors corresponding to the boundary flow 
(—)&(+), as follows from 



#(-o)fe(o+) 



^(1 + ^5) 



(5.29) 
-0)&(0+) -> 

(5.30) 



9{-)H+) 2 

A plot of In g in ( 5.23 ) as a function of the boundary scaling parameter 9b defined in ( |5.25|) 
with finite n P] for various values of r is given in Fig. |3|. Observe that for r << 1, the correct 
conformal boundary entropy is reproduced. As r increases, one can see that the entropy 
deviates from the conformal field theory value. Indeed, for r — > 1, the entropy approaches 
zero, as expected for a massive theory. 



10 



The horizontal axis is rescaled in such a way that the range is mapped to (0 , 1). 
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boundary entropy vs. boundary scaling parameter 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 3: 



One might wonder how there can be a flow to the boundary condition (—)&(+) in the 
N = even - Sector II, for which the boundary "spins" are fixed to (see ( [4.2[ )). Our 
explanation is that there are boundary bound states with spins ±1, corresponding to the 
pole at = O = 9 B - f in the CDD factor. (See Figure |.) 



5.3.2 The boundary flow (-0)&(0+) -> 2(0) 
Let us now consider instead the scaling limit 



71 

£ = i(0B-ln-), n^0, (5.31) 



with 6b finite. Taking again the sign — in the limit ( gjg) , the factor P CDD (9 1 i) becomes 



real, and so the corresponding kernel k^{0) vanishes. However, 

, isinh(^-^) 

Q(8,0^—-yfj (5-32) 

sinh(^ + f ) 

and therefore 

^(0)->$(0-0 B ). (5.33) 
The result ( |5.23| ) for the boundary entropy now implies 



\ ng = —l d6 $(0 - B ) L 2 (6) . (5.34) 
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Figure 4: Boundary bound state pole. 



Using the results L 2 (— oo) = In (|(3 + \/5)), 1/2(00 ) = In 2, we obtain Q 



This is the ratio of g factors corresponding to the boundary flow (— 0)&(0+) 

fl , (-0)&(0+) 



#2(0) 



(5.35) 
2(0), since 
(5.36) 



5.3.3 Nonsupersymmetric flows 

The analysis presented so far in Sees. f| and |5| has been restricted to the NS case of the 



TIM, for which the boundary S matrix is given by ( 2.21 ), ( 2.22 ). However, the results for 
the cases without supersymmetry can now be obtained with no additional effort. 

For defmiteness, let us now consider the nonsupersymmetric boundary S matrix fl2.12| )- 
( |2.18 ). Q The corresponding inversion identity is again given by ( |4.6j )-( p~8| ), except the 



11 This flow does not occur for N = even - Sector I, since for this sector there is no £2(0) contribution to 



the boundary entropy, as can be seen from (5.21). Our understanding of this fact is as follows: By definition 



(4.2), this sector has boundary "spins" ±1. Moreover, in the scaling limit (5.31), there cannot be a boundary 
bound state with spin 0, since the CDD factor does not have a corresponding pole. That is, the process 
represented by Figure^ with the spins ±1 and interchanged does not occur. Hence, there cannot be a flow 
to the boundary condition (0) in this sector. 



12 For the other nonsupersymmetric case (2.19), the results are exactly parallel, with the spins —1 inter- 
changed with +1. 
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sectors are now given by 

N = even — Sector / : a± = a[ — a^+i = a' N+1 = — 1 



N = even — Sector II : a x = a[ = a^+i = o! N+1 = 



iV = odd — Sector / : a\ = a\ = — 1 , O/v+i = a N+i = 



N+l — u ) 

/ 

N+l — w j 



N = odd — Sector 77 : a\ — a x — , cln+i = o-'n+i = ~ 1 • (5.37) 
That is, the sectors are restrictions of those in the NS case (f4.2p. In particular, the iV = 



even - Sector II is identical to the one for the NS case. Hence, the TBA equations and 
boundary entropy are the same as before ( |5.19|) , ( |5.23|) . Moreover, the two massless scaling 



limits give the same results (|5.28| ), ( |5.34| ). However, the interpretation of these scaling limits 



is different from the interpretation in the NS case: the first scaling limit now corresponds 
to the boundary flow (—0) — > (— ), while the second scaling limit now corresponds to the 
boundary flow (—0) — > (0). That both interpretations are possible is due to the coincidence 
in the ratio of g factors 0, 

0(-o)fe(o+) = #J-o^ ff(-o)fe(o+) = £j-o^ ^ 3g ^ 

£(-)&(+) 9{-) ' 5-2(0) 9(o) 

The TBA results (|5.28| ), (|5.34j ) for these flows coincide with those obtained in [31| on the 



basis of an analogy with the Kondo problem. 



6 R case 



We now consider the R case of the TIM, for which the boundary S matrix is given by 
fl2.24j )-( p.26| ). Remarkably, the results are closely related (in fact, dual) to those for the NS 
case. Indeed, let us define the four sectors as before (|4.2|) . The relation between the full and 
reduced transfer matrices is again given by (^4.3|) , except w^ a \9) is now given by 



N 



W 



(a) 



{ N(i7r-e,£ + ) N(e,£j) 

R(i7r-9,Z+) R(6,£J) 
{ e 2 ^ e R{m-9^ + ) N(9,£..) 



The inversion identity is again given by ( |4.6| ), (|4.7|), except f± (9) is now given by 

^(cosh^ ± cos £_) (cosh 9 ± cos £ + ) 



J± 



(d) 



4 

i(cosh 9 ± cos£ H 
2 (cosh 9 ± cos £_ 



(6.1) 



(6.2) 
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The periodicity and crossing properties of the reduced transfer matrix are the same as before 
(PI), ( Ei~l~0D ■ In contrast to the NS case, the reduced transfer matrix now becomes an anti- 
diagonal (rather than diagonal) matrix for 9 — > oo. Nevertheless, the asymptotic values of 
the eigenvalues are again given by ( |4.17|) , except z^ a > is now given by 



(0) 



± - 



±2 (-if 



±[-f 



-+1 



±2 



ie 
4 



JV+l 
2 



JV+l 
2 



(6.3) 



A suitable Ansatz for the eigenvalues is again given by ( ELT8D , except c (q) and S a) are now 
given by 



f ±1 

±2 
±1 
±2 



IV 



+ 1 



N 
2 

jV+1 
2 

N+l 



(6.4) 



respectively. The Bethe Ansatz equations are therefore again given by ( |4.20| ), with the new 
f±\0) given in (|6.2|) . Comparing with the old f± (9) given in (|4.8|) , we conclude that the 



Bethe Ansatz equations for the R case exactly coincide with those for the NS case, except 
the sectors I and II are interchanged (for both N = even and N = odd)! We remark that 
the eigenvalues do not depend on the parameter r which appears in the boundary S matrix. 

It is now straightforward to repeat the TBA analysis. For N = even - Sector I, we obtain 
the same constraint equations (|5.5|), ( |5.13| ), and therefore the same TBA equations fl5.19| ) 
and boundary entropy ( |5.23| ). The result ( |5.28|) for the first massless scaling limit can now 
be interpreted as the boundary flow (d) — > (0), since 

M = 7 f l(i + V5), (6 . 5) 
9(0) 1 

as follows from ( |2.8| ). Similarly, the result ( |5.34j ) for the second massless scaling limit can 
now be interpreted as the boundary flow (d) — > (—)&(+), since 

- 2 1 / 

3 + ^. (6.6) 



9(d) 



9(-)k(+) 



V2 



7 Conclusion 



We have achieved the principal goals set out in the Introduction: 
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We have provided support for the proposed TIM boundary S matrices M, [| by show- 
ing that the corresponding boundary entropies ( |5.23|) , (|5.28|) , ( {5.34T ) are consistent with 
boundary flows (both supersymmetric ( |5.30| ), ( |5.36| ), (|6.5| ), ( |6.6| ) and nonsupersymmet- 
ric ( |5.38| )) which were expected on other grounds 
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We have developed in Appendix || analytical tools for RSOS models with boundary, 
which we have used to derive exact inversion identities for the TIM. (See ( |4.6| )-( P~E[ ) 
and ( |6.2|) for the supersymmetric cases, and ( |5.37|) for the nonsupersymmetric case.) 



Our TBA descriptions of boundary flows have been derived directly from the TIM 
scattering theory. The fact that we have reproduced the TBA description of the non- 
supersymmetric flows given by Lesage et al. |31| provides support for their approach 
based on an analogy with the Kondo problem. The TBA descriptions of the super- 
symmetric boundary flows are new. 



It would be interesting to see if the approach presented here can also be used to investigate 
massless flow in the bulk |H] . Moreover, we expect that it should be possible to generalize 
this approach to more complicated models, such as the RSOS n models with n > 3 BI| , 
and coset models JIT], 
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A Properties of S matrices 

We collect here some important properties which are satisfied by the TIM bulk and boundary 
S matrices. 



A.l Bulk S matrix 

The bulk S matrix ( |2.3|) - fl2.6|) has the following symmetries in its indices 

S c a d b (0) = S c b d a (9) = S d a c b (9)- (A.l) 
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It also satisfies the crossing relation 

S c a d b (9)=S a c b d (iir- 

and the unitarity relation 



E 5 a6(0) S d a t{-9) = 5 cc ,A ac A bc . 



where A a b is the so-called adjacency matrix (see e.g. p9| ) 

A a b = $a,b-l + o ajb+ i . 
Moreover, it satisfies the Yang-Baxter (star-triangle) equation 

E s « oiOi - o 2 ) si *{0i - e 3 ) si i{9 2 - *s) 

9 

9 

Finally, we note that the bulk S matrix at zero rapidity is given by 

S c a d b (0) = 5 c4 A ac A bd . 



(A.2) 
(A.3) 

(A.4) 



(A.5) 



(A.6) 



A.2 Boundary S matrix 

The non-supersymmetric boundary S matrix fl2.12j )-( [2.i8| ) obeys the unitarity relation 

Y^Ra C (6)Ra d (-6) = 6 b4 A ab B d , (A.7) 



where here B d equals 1 if d is an allowed state of the boundary and equals zero otherwise; 
hence, it is given by 



Bd = 8d,-i + $d,o ■ 
This S matrix also obeys the boundary crossing-unitarity relation M 



Rb 



£«2^<(f + 



(A.8) 



(A.9) 



as well as the boundary Yang-Baxter equation [BL [13], f2£J, [K], [39 



E S » ^ - 2)R9 f (6,) S 9 a fa + 9 2 )Rd \ (9 2 ) 



f,9 



(A.10) 



= E^ b (e 2 ) s b a f g (9i + e 2 ) Rf (e,) si d e {e x -e 2 ). 

; C 9 

f ,9 

The supersymmetric boundary S matrices described in Sec. |2.2.2| obey the unitarity condi- 
tion ( |Q|) with B d = 1, and also flPD , flOQ|) . 
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B Fusion procedure for RSOS models with boundary 



The fusion procedure was developed for bulk vertex models in |25 
to bulk RSOS models in 



26f, and was adapted 



The fusion procedure was extended to vertex models with 
boundary in |?7|], but this work has been adapted only in part to the RSOS case |29], |30 



In particular, the useful notions of projectors and quantum determinants have not been 
explicitly implemented in [^, |3(| . For this reason, and also to make this paper self-contained, 
we give here a brief summary of the fusion procedure for RSOS models with boundary, 
and provide the derivation of the TIM inversion identity. However, our treatment is not 
completely general. In particular, to avoid complications which are not necessary for the 
TIM, we restrict to S matrices with the symmetries (|A.1|) . 

We remind the reader that a bar over a quantity (e.g., S) denotes that it is "reduced," 
and a tilde over a quantity (e.g., S) denotes that it is "fused." 



B.l Projectors 



We shall carry out the fusion procedure by exploiting the fact that the reduced bulk S 

c d 
a b 



matrix degenerates into the projector P c d h for some value of the rapidity, which for the 



TIM is 9 = -in, 

S: d b {-m) = V2P-l d b . (B.l) 
For the TIM, P~ c a d b has the matrix elements 

P-li. = *«,-*>, (B-2) 
where as usual a ,a' £ { — 1 , +1}. We define the projector P + c a f by 

P +C ai = Kt-P- C ai, (B.3) 

where l c a f is the "adjacency- inclusive" identity matrix, 

F a d b = 5 c>d A ac A bd . (B.4) 

13 Note that we work here with the reduced matrix f(0) rather than the full matrix t{&). There are 
good reasons for so doing: (i) as explained in Sec. it is the reduced transfer matrix for which we require 
an inversion identity; and (ii) the full S matrix is singular at 9 = —in. 
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k 


c 


d 



62 



= c 3 



Figure 5: Fused bulk S matrix S b k 



a d 



For the TIM, P + c a d h has the matrix elements 



P+ O" CT 

r 



p+ 
r a a' 



^OW' + $a,-<T> 



1 

2 ' 



<5o-cr' • 



These projectors have the important properties 



6 



Ep- c c' p- c' 

r a b r a 

d 

Ep- c d p+ d 

r a b r a 



p- c a 

r a b ' 



0. 



2> 



;b.5) 



+ c c' p+ c' d 
a b a b 



c d 



[B.6) 



B.2 Fused bulk S matrices 



We derive a bulk "fusion identity" from a degeneration of the Yang-Baxter equation. That 
is, in (|A.5| ) we set 9\ = 9, 9 2 = 9 + in, 9 3 = 0, use the degeneration result (|B.1|) , and contract 
on the right of both sides with the projector P + to obtain 



(B.7) 



f,9 

This identity can be used to show that the "fused" S matrix (which can be read off from 
by replacing the projector P~ with P + , and which is represented by Figure [5p, 



Bb a k d w = E p+ -^)^+m)p 

f,9 



+ f k 
j d 



(B.8) 



satisfies the generalized Yang-Baxter equation 

- b f 



J2Sk 9 (9 1 -9 2 )S 9 i (9,- 9 3 ) S[l(9 2 - 9 3 

' ' a c J d 

f,9 



f,9 



S b 3 
k g 

a f 



% d(P* -dz)Sk -g (9, - 9 3 ) S > \ (9, - 9 2 ) 



~k. f m 
g i 

j d 



(B.9) 
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~/b k d 

Figure 6: Fused bulk S matrix S a c • (6\ — 9 2 ) 



For the TIM, the nonzero matrix elements of S are given by 

~ cr -Oct I 

S * (9) = S^o (9) = — cosh - 

v2 2 



cr 



cr 



From a second degeneration of the Yang-Baxter equation ( |A.5| ) with 9 2 — 9 3 
obtain a second fused S matrix (see Figure §) 

/>9 



which obeys 



f,B 



% f 3 



J2 S \{ (9 2 - 9 3 ) S b a i f (9 t - 9 3 ) S f s l a d m {9i - 9 2 ) 



f,9 



(B.10) 



— iir, we 



(b.ii; 



(B.12) 



For the TIM, the nonzero matrix elements of S are given by 

~/0 ~la cr' cr 1 

s * cr' A°) = S (9) = -j= cosh - . 



(B.13) 



B.3 Fused boundary S matrix 



_ k 
d 

f 



;b.i4) 



Following [27], we obtain a boundary fusion identity from the degeneration of the boundary 
Yang-Baxter equation ( |A.10|) with 9\ — 92 = —in, 

d,f,g 

This identity can be used to show that the "fused" R matrix (which can be read off from 
( |B.14| ) by replacing the projector P~ with P + , and which is represented by Figure [?]), 

~ R \ k W= E P+ -^ ff/ (0)S 9 a d f (29 + i7r)Rd k (9 + i n)P +d a i (B.15) 

c *— ' c 1 f 

d,f,g 
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- j k 

Figure 7: Fused boundary S matrix R a {9), 

b c 



satisfies the generalized boundary Yang-Baxter equation 

f ,9,h,i 



x — , -la he f - g e - j k 

E s dbg (e 1 -e 2 )R9 (e 1 )Shi(e 1 + e 2 )Re (e 2 ) 



E~ i g -la h g _ k ~ f e 

Ra (9 2 )S dl f (9 1 + 9 2 + m) Rf (^) S h j (9 1 - 9 2 - in) . (B.16) 
be 1 g d k 

f ,9,h,i 

The shifts in the arguments of the fused bulk S matrices on the RHS should be noted. As 
an example, for the NS case of the TIM, the nonzero matrix elements of R are given by 

Ra° ° (0 f) = v^cosh-, ^o" ° (9,n = ^ (cosh + cos cosh-. (B.17) 
o a v s/ 2 a' o v ' 2y/2 2 

B.4 Fused transfer matrix 

Before attempting to construct the fused transfer matrix, it is instructive to first review the 
construction of the fundamental transfer matrix. To this end, we set 



u,± i • • •'J'iV -|- 1 

R + ,«? (0-i7T,£ + )T , a" (9\9 1 ,...,e N ), (B.18) 

•• 1 - ., a l °1— °iV+l 

a'{ 

where T is defined by 

r ai ;'a( +1 {e\e u ... f e N ) = £ { f t'-ttl W > ■ ■ ■ > ^) £ <r (0 , £_) 

a i- a iv+i „ „ L 1 JV+1 a iV+l 

x r m , . . . ,e N )} , (B.19) 
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and the monodromy matrices T and T are given by 

f m,-, N ) = < J (9-9,)... S% X + + l (6 - o N/ 

T XZ + + \ m,-.., 9 N ) = (* + '")••• S% 2 + ( B - 2 °) 

The boundary matrix R in ( |B.19| ) is assumed to obey the boundary Yang-Baxter equation 
( |A.10| ), which implies that T obeys 

E_ _ 6i...bjv+i _ _ ai...a N+ i 

Sl{\ (0i ~ 02) T * (6,) S C J * (0 a + 9 2 ) T * (6> 2 ) 
i i a 1 ...a JV+1 l bi...6jv+l 

ci ,6i ,Ojv+i 

= e f TT ^ + ^ rT +1 ^ - ^ • (R21) 

a l--- a jV+l 1 01---0JV+1 1 

ci ,6i ,0jv+1 



However, the matrix R + in ( |B.18|) is not yet specified. Indeed, following Sklyanin f24 |, the 



requirement that the transfer matrix obey the commutativity property ( |3.14| ) determines the 
relation which R + should satisfy. In this way, we find (using also the properties (|A.1|) - ([X73|) ) 
that 

R +C a (9,{ + ) = Ra C (B.22) 

b b 

where R a (9,£) obeys ( |A.10| ). The result ( |B.18| ), ( |B.22| ) coincides with the expression 
for the (reduced) fundamental open-chain transfer matrix . 

We follow a similar strategy to construct the fused transfer matrix i. We set 

fl *)= E R , >" (fl-^.f^W (9\9 U ...,9 N ), (B.23) 

where T is defined by 

~ ai...a N+ i _ — ^ ( ~ a l ...a N+1 

T a» ,6" ,c? (0|0 X , . . . , ^) = V \ T K K +1 (9\0! , . . . , 9 N ) 

a 2 ,a JV + 1 , 

6" c" 
JV+1 ' JV+1 

~ c jv+i a iv+i X ai...ajv+i 

x ^ (9,£_)T 4,c^ +1 (e + iir^,... ,9 N )\, (B.24) 

fe jV+l a JV+l a l--- a jv + i ■> 

where the fused monodromy matrices T and T are given by 

~ Ol...Ojv+l X — -\ - a 'l a 2 - a JV a iV+l 

T &! ,bN+i (9^ , . . . , 9 N ) = ) S hb 2 (9-9t)...S b N b N+1 (6 ' 

a[...a' N+1 ai a' 2 a N a' N+1 

O2,— ,0 N 



'N ) 



X ai...apf+i x — — / A'.: I i bjv o-N+l -laL b\ 0,2 

T 6; (^,...,61^)= «* (* + • • " 5 < b 2 ai (d + Ox) (B.25) 



62 ,— i b IV 
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We determine the relation obeyed by R from the requirement that the fused transfer 
matrix ( |B.23| ) commute with the fundamental transfer matrix ( |B.18j ), ( |B.22 ), 



t(9\0i , • • • , 9 N ) , t(9'\9i , . . • , 9 N ) 
With the help of the relation obeyed by X 



0. 



E 

g ,h ,i , 



~/a h ai - a-, ...a N ,, ~ g e ~ bi...bjv+i 

S d (9 1 -9 2 )T 9 + (9,) S h i (9, + 9 2 ) T f e,ij (9 2 ) 



E 

g ,h,i, 



(B.26) 



, , -la ha', - b\...b N+ i ~ / e 

T a,bf (9 2 ) S di f \9 l + 9 2 + t7r) T f (9,) S h j (9 1 - 9 2 - iir) 



(B.27) 



we obtain the following equation for R 

E~ g k _ d -la f d ~ + e c 

Sfi (9 2 -9 1 )R + B (e 1 -m,Z + )S jbg (-9 1 -9 2 + 2m)R a(0 2 -nr,£ + ) 

7 m m d b 

b,d,f,g 



E~A-db ~ g k _ e ~/o f e 



b,d,f,9 



3 d 



(B.28) 



That is, this relation guarantees the commutativity ([B.26|). This relation is satisfied by 



~ + k j 
R 

c b 



£ + )= E^-^ (^-d^ + )S 9 f d a (-29 + in)Rd k (-9,£ + ) P +d k {. 

z — ' c J f 



dj ,g 



(B.29) 



For the NS case of the TIM, the nonzero matrix elements of R are given by 



~ + a 

R 

a 



U <T 



\/2 cosh - , 

1 $ 
— — (cosh 9 + cos £+) cosh - 
2v/2 +; 2 



(B.30) 



To summarize, the fused transfer matrix i is given by (|B.23|) -( [B.25|) , where the fused 
matrices S, S , R and R are given by ( B.8|) , (B.ll| ), (B.15| ) and (|B.29|) , respectively. For 
the NS case of the TIM, the nonzero matrix elements of the fused transfer matrix are as 
follows: for N = even, 



t 



pi CT 2 ... a N Q 

a l +1 (9\9 1 ,...,9 N ) = 2 cosh 2 - (J] cosh), 



CTl <T2 ... 



(B.31) 



~0 (Tl CT2-.CT JV 
^0 (Tl CT2 



■ CT JV 1 2® TV 

^ (9\9i , . . . , 9n) = - (cosh 9 + cos £+) (cosh 6* + cos £_) cosh -(J^cosh) 
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and for N = odd, 

_o-i a 2 ...a N+1 Q 

o «...<r^ Q {B\e x ,... ,9 N ) = (coshfl + cos £_) cosh - (jjcosh), 

to* 1 o...o<,l^(0\0i,---M = (coshtf + cos £ + ) cosh 2 - (Q cosh), (B.32) 
where (Y[ cosh) denotes 

N 1 1 

(Jjcosh) = Y[cosh(-(8-6 j ))cosh(-(6 + 6 j )). (B.33) 

3=1 

As also discussed in Sec. 0L for a given a transfer matrix (either fundamental t™, 1 '" f +1 or 

i_i a 1 ...a i y. +1 

fused £„/ _/ ), it is convenient to define the following four "sectors": 

AT = even — Sector / : ai , a' x , a N+ i , a^ +1 G { — 1 , +1} , 

N = even — Sector 77 : a\ — a'± = a^+i = a' N+1 = , 

N = odd — Sector/ : a\ , a[ G { — 1 , +1} , ajv+i = o,' N+1 = , 

N = odd — Sector// : a\ = a[ = , ajv+i , a' N+1 G { — 1 , +1} . (B.34) 

The results (|B.31| ), (|B.32|) show that, within each sector, the fused transfer matrix is pro- 
portional to the adjacency- inclusive identity matrix, 

f a \e\e 1 ,...,e N ) = g< a \e)ito, (B.35) 

where a runs over the four sectors ( p.34j) , and g^ a \0) is given by 



2 

\ (cosh 6 + cos £ + ) (cosh 6 + cos £_ 



0.36) 

(cosh^ + cos£_) 



g^\e) = cosh 2 °- f[ cosh(i(0 - 9 3 )) cosh(^(0 + 6 3 )) x < 

3=1 

(cosh 6* + cos£ + ) 

respectively. This is a nontrivial property of the TIM. The supersymmetric sinh-Gordon 



model enjoys [21] a similar property. 



B.5 Fusion formula and quantum determinants 

We now derive the important "fusion formula," from which the TIM inversion identity is 
obtained. To this end, we first note that T ( |B.24| ) can be expressed as the fusion of the 
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corresponding fundamental quantities T ( |B.19| ), 

~ ai...a N+1 _ — ^ f , ui f _ bi...bf 



~ a\...a N+ i _ — ^ r , f _ t)i...0jv+i _r 

r < m A [B\ Bx , . . . , Bn ) = £ {p + * { r a , . . . , ^) s§ g (20 + 

a 1 ...a JV+1 ^ 11 a i--- a jv+i 1 



Z7T 



i>i >f>jv+i : 

/l .91 



_ ax ...ajv-i-i r „ 1 

x T 9 i (0+i7r|0i,...,^)P + ^5,[. (B.37) 
We next observe that the reduced bulk S matrix obeys 

^TS^V - 29) S: b d ,{iir + 26) = 6 d>d/ Q{e)A a d A bd , (B.38) 

c 

where, for the TIM, the scalar factor Cd{6) is given by 

£ o (0) = icoshfl, C ±1 (0) = 2 cosh 0. (B.39) 



Using also Eqs. ( |B.29| ) and ( |B.21| ), we obtain the desired fusion formula Q 



•j" n" 

1 ' U JV 



where the quantum determinant A(0) of the transfer matrix is given by 



/l ■■■■ ,/iV+l 
b ,c ,d ,j ,k 



x )); ^rif + + : *(r (*))} h'-'-'-fN+i } > (b- 41 ) 

where the quantum determinants of the monodromy matrices are defined by 

5(p(0))° a , 1 --- a f +1 = V p- c f^-^W) f^- h aN+1 (e + in) P~ bN+1 \ , 

V V lib a' v ..a' N+1 / j a\ a 1 a 1 ...a N+1 \ I bi...b N+1 \ / a N+1 a' N+1 ' 

61 ,fejv + i 

6(T(6)) e h a rJ +1 = V P~ c , bl t')'"^ +1 (6) f ai - aN+1 (0 + m ) p- \ 6 ^+i (9.42) 

V V lib a' v ..a' N+1 / ■> a\ ax a i-- a N +i K ' bx—bN+x v / a jv+i a iv+i ' 

and the quantum determinants of the boundary matrices are defined by 
5(R(9,^))a jk = ^P-^P/ {9)Slpe + m)R<i k {e + i<K)p- d a i, (B.43) 

b c ^— ' c J f 

d,f,g 

5(R+(6,Z + )) kj a = ^ P ~c^Ra\^-0^ + )Sj d a (-2e + l n)R d k (-9,^)p- d k i. 

d,f,g 

1 We save writing by suppressing the dependence of the transfer matrix, etc. on the inhomogeneity 
parameters 61 ,9n- 
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We now proceed to evaluate the quantum determinants for the TIM. With the help of 
the identity 

£ P ~- S c 9 1(0) SI }{9 + in) P- \ k d = iy/2 sinh °- P~ \{ P~ ° * , (B.44) 
f,g 

we find that the quantum determinants of the monodromy matrices are given by 
8(f(9))t = (iy/2 sinh -) N P~ c a ? P~ a[ T . . . P~ a ' N b , , 

v V lib a' v ..a' N+1 \ 2 a '2 a-2 a-N+i a N+1 ' 

5(T(9)) C , a , 1 -T +1 = (iV2 sinh -) N P~ c 1 P" ai ^ . . . P~ ajv * . (B.45) 

v v //b a' v ..a' N+1 \ 2' a i °l a 2 a 2 a iv+i V ' 

Moreover, for the NS case, the quantum determinants of the boundary matrices have the 
following nonzero matrix elements 

8(R(9,£-)) -° ° ° = zv^sinh^, 

a 2 

5{R{9, £_)) = ± -L= (cos - cosh 0) sinh Q -, a = ±a' ; (B.46) 

and 

5CR + (0,f+)) CT °- CT = -iV^rnnh^, 

CT 2 

(5( J R + (0,e + )) o °J = ±-^(cosh£-cos£+)sinh^, a = ±a' . (B.47) 

We conclude that the quantum determinant A(#) for the NS case of the TIM is given by 

A (a) (#) = h ia) {6) , (B.48) 
where a runs over the four sectors ( p.34j) , and h^ a \8) is given by 

2 

9 N 



h^(9) = sinh 2 - J] sinh(^ - 6 3 )) sinh(^ + 0,-)) 



i 



\ (cosh — cos £ + ) (cosh 9 — cos£_) 



V 

2 '2' ' '2' J (cosh — cos £ 

(cosh0 — cos^_) 



respectively. Substituting this result, together with the result ( p.35|) , (|B.36|) into the fusion 
formula ( B.40|) , we finally arrive at the inversion identity (|4.6|) - (|4.8|) . 
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